High Density Neutron Star Equation of State from 4U 1636-53 Observations 
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A bound on the compactness of the neutron star in the low mass x-ray binary 4U 1636-53 is 
used to estimate the equation of state of neutron star matter at high density. Observations of 
580 Hz oscillations during the rising phase of x-ray bursts from this system appear to be due to 
two antipodal hot spots on the surface of an accreting neutron star rotating at 290 Hz, implying 
the compactness of the neutron star is less than 0.163 at the 90% confidence level. The equation 
of state of high density neutron star matter estimated from this compactness limit is significantly 
stiffer than extrapolations to high density of equations of state determined by fits of experimental 
nucleon-nucleon scattering data and properties of light nuclei to two- and three-body interaction 
potentials. 

PACS numbers: 26.60.-l-c, 97.60.Jd, 21.65.-ff 



I. INTRODUCTION 

Oscillations in x-ray brightness during bursts have 
been observed from ten low mass x-ray binaries with 
the Rossi X-ray Timing Explorer 0]. These oscillations 
appear to be due to rotational modulation of one or a 
pair of antipodal hotspots on the surface of the accreting 
neutron star With this model the amplitude of the 
oscillations is determined by the mass-to-radius ratio of 
the neutron star, the compactness M/R. A more com- 
pact neutron star results in greater deflection of the x-ray 
photons emanating from a hot spot, making the hot spot 
visible to an observer for a larger portion of the rotation 
period and so reducing the amplitude of the brightness 
oscillations. The bounds on the compactness of the neu- 
tron star are particularly limiting when two antipodal 
hot spots are present [Q. 

The source 4U 1636-53 appears to be the best candi- 
date for providing a stringent limit on the compactness of 
neutron stars from x-ray burst oscillations |4| . A 290 Hz 
subharmonic of the stronger 580 Hz brightness oscillation 
has been observed in five bursts from this source j5j , but 
the subharmonic has not been confirmed in subsequent 
bursts . The existence of the subharmonic would sug- 
gest the spin frequency of the neutron star is 290 Hz, 
and the presence of two antipodal hot spots is responsi- 
ble for the 580 Hz modulation. A pair of high frequency 
quasiperiodic oscillations separated by 251 Hz has also 
been reported for this source and interpreted in a beat 
frequency model to be a consequence of a spin frequency 
of 290 Hz with two hot spots on the neutron star ^ . The 
x-ray oscillations during the rising phase of bursts from 
4U 1636-53 have been recently modeled as due to a pair of 
circular antipodal hot spots, each having an angular size 
that grows linearly in time |4j. The two hot spot model 
constrains the neutron star compactness to M/R < 0.163 
at 90% confidence (using G — c— 1), requiring the radius 
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of a 1.4 Mq neutron star to be greater than 12.7 km. The 
99% confidence level compactness limit, M/R < 0.183, 
implies i? > 11.3 km when M = 1.4Af0. 

During a burst event in low mass x-ray binaries the os- 
cillation frequency is usually observed to evolve asymp- 
totically to a limit that probably corresponds to the spin 
rate of the neutron star. Measurement of the asymp- 
totic frequency at different orbital phases of the binary 
could reveal a Doppler modulation that would allow for 
a determination of the masses of the neutron star and 
the companion. A mass measurement together with the 
compactness limit would tightly constrain the radius of 
the star. A Doppler modulation has not been observed 
for 4U 1636-53 in an analysis of 26 burst oscillations [Q. 
The absence of the modulation likely limits the neutron 
star mass in this system to M < I.6M0 Q. 

The purpose of this paper is to use the 90% confi- 
dence level compactness bound on the 4U 1636-53 neu- 
tron star to estimate the high density equation of state 
of neutron star matter. Ideally one would have both a 
measured mass and radius for a number of neutron stars 
over a wide range of masses up to the maximum mass 
value. The observed mass-radius relationship could then 
be used to uniquely calculate the equation of state by in- 
verting the Oppenhemier-Volkoff equation, the relativis- 
tic stellar structure solution to Einstein's equations 
Even if the mass and radius are known for only a sin- 
gle neutron star, the high density equation of state can 
still be accurately estimated up to the central density 
yielding the mass of this star. At low density the equa- 
tion of state can be determined by fitting experimental 
nucleon-nucleon scattering data and properties of light 
nuclei to two-body and three-body interaction potentials 
1^, . Confidence in these low density equations of state 
is high near normal nuclear density (energy density =153 
MeV/fm"^, baryon density = 0.16 fm"'^) because nuclear 
physics experiments allow for verification in the low den- 
sity regime. Estimates of the maximum possible mass 
of neutron stars based on these equations of state have 
regarded them as experimentally verified up to twice nor- 
mal nuclear density ^1, 12 1. Assuming one has values of 
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both the mass and radius for a single neutron star, the 
equation of state in the density range between the experi- 
mentally verified low density regime and the high central 
density of this star can be accurately estimated to have a 
polytrope form provided the matter does not undergo a 
phase transition in this density range |^ . The polytrope 
high density equation of state is chosen to match the en- 
ergy density and pressure of the low density equation of 
state at the upper limit of applicability of the low den- 
sity equation, and the adiabatic index of the polytrope 
is chosen to yield a stellar radius corresponding to the 
observed radius of the star when the mass of the model 
star equals the observed mass. The resulting polytrope 
is an estimate of the equation of state of neutron star 
matter in the high density regime. 

In the case of the neutron star in 4U 1636-53 we do not 
have measurements of both its mass and radius, but only 
an upper bound on the compactness and an upper limit 
on the mass. Nevertheless the high density equation of 
state is still tightly constrained by the compactness limit 
provided it is assumed the neutron star mass is at least 
1.4M0, a plausible assumption for an old accreting sys- 
tem like 4U 1636-53. In the next section the high density 
equation of state is estimated from the 90% confidence 
level compactness bound on the neutron star in 4U 1636- 
53 assuming the neutron star mass is at least 1.4Mq. 
The main result of this paper is that the A18 + Sv + 
UIX* equation of state of Ref. jl^ (hereafter the APR 
equation of state), a recently developed equation of state 
representative of the class of equations of state derived 
from fits of experimental nuclear physics data to inter- 
action potentials, is too soft at high density to yield a 
star with the observed low compactness of the neutron 
star in 4U 1636-53. Further if the APR equation of state 
is used up to twice normal nuclear energy density, and 
any polytrope equation of state having a sound speed 
less than the speed of light is used above twice normal 
nuclear density, the model star compactness still exceeds 
the compactness of the 4U 1636-53 neutron star. Using 
the APR equation of state up to just a value in the range 
1.55-1.64 times normal nuclear energy density and a stiff 
polytrope with an adiabatic index rof5.5^r<7at 
high density can produce a star with the low compact- 
ness of 4U 1636-53 without the speed of sound waves in 
the neutron star interior exceeding the speed of light. 



II. ESTIMATION OF THE HIGH DENSITY 
EQUATION OF STATE 

The equation of state of high density neutron star mat- 
ter is estimated in this section from the compactness 
bound on the neutron star in 4U 1636-53. The equation 
of state at low density is assumed to be known for energy 
densities less than a transition density p. Here the APR 
equation of state is used for the low density portions of 
the stellar interior, from normal nuclear density p„„ up 
to p. The BPS equation of state is used for densi- 



ties below pnn to model the matter outward through the 
neutron star crust to the surface. Recent work on deter- 
mining the maximum possible mass of neutron stars has 
assumed the APR equation is accurate up to twice nor- 
mal nuclear energy density 1 1 1 , 12] . For densities greater 



than p an equation of state of polytrope form is used in 
this paper: 
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The polytrope form can be expected to be a good approx- 
imation to the actual high density equation of state of 
the 4U 1636-53 neutron star provided the matter has not 
undergone a phase change to a softer condensed phase, 
such as a kaon condensate or a deconfined quark phase, 
at some density above p The compactness limit on 
the 4U 1636-53 neutron star requires the stellar radius to 
be large, so it is unlikely this star contains a condensed 
phase. 

The algorithim for determining the high density equa- 
tion of state begins by choosing the transition energy 
density p where one switches between the high and low 
density equations of state. The pressure p at the transi- 
tion density is computed from the low density equation of 
state and the transition energy density value: p — p (p). 
The adiabatic index P in the polytrope high density equa- 
tion of state Eq. (|l|) is treated here as a free parameter 
that is varied from small values (soft equations of state) 
to progressively larger values (stiffer equations of state) 
to find the set of model stars satisfying the 4U 1636-53 
neutron star compactness bound. 

The interior structure of each model star is assumed 
to be described accurately by the Oppenheimer-Volkoff 
solution, the solution to Einstein's equations for the in- 
terior of a nonrotating star in general relativity: 



dp {p + p) [m{r) + iirr^p] 
dr r[r — 2m{r)] 



where 
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Rotational corrections to the Oppenheimer-Volkoff solu- 
tion are small except for very rapidly rotating stars 
and will be ignored here. Families of model stars, one 
family for each (p, F) pair, are determined by integrat- 
ing the Oppenheimer-Volkoff solution Eq. (g) for a set 
of central energv density values. The polytrope equation 
of state Eq. (H) is used for densities exceeding p. The 
APR equation of state is used for densities below p and 
above normal nuclear density. The BPS equation of state 
is used below p„„. 

Figure |^ shows how the neutron star mass depends on 
the radius for the case where p = 2p„„. Several families 
of model stars are represented differing by the value of 
the adiabatic index of the high density polytrope. Also 



a 



2.0 



1.5 



1.0 



0.5 



r = 3 




APR \\ 


\\ 










'/^M/R = ai83 / 


V / 




1 / 


//' M/R = 0.163 . 


r = 2 1 / 1 




1 





10 



11 12 13 

Radius (km) 



14 



FIG. 1: Neutron star mass vs. radius for the APR equa- 
tion of state used up to twice normal nuclear energy density. 
The solid curves are for a polytrope high density equation of 
state, and the dashed curve is for the APR equation of state 
used at all densities. None of these models satisfy the 90% 
confidence level 4U 1636-53 neutron star compactness limit, 
M/R < 0.163. 



shown is the neutron star mass-radius curve resulting 
when the APR equation of state is used at all densities. 
Lines representing the 90% and 99% confidence level 4U 
1636-53 neutron star compactness bounds delineate the 
allowed regions of the mass-radius parameter space. The 
allowed regions are below these compactness limit lines 
and above the IAMq line. The compactness of all these 
model stars exceeds the 90% confidence level 4U 1636-53 
neutron star compactness bound, M/R < 0.163, even for 
very stiff polytropes. Models with F > 7 are not physi- 
cally realistic because the equation of state predicts the 
speed of sound waves. 
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FIG. 2: Neutron star mass vs. radius for the APR equation 
of state used up to 1.64 times normal nuclear energy density 
and a polytrope equation of state used for higher densities. 
Models with F ^7 satisfy the 90% confidence level 4U 1636- 
53 neutron star compactness limit, M/R < 0.163, without 
the sound speed exceeding light speed. 



TABLE I: Equation of state of neutron star matter using 
the APR equation of state up to 1.64p„„ and a F = 7.023 
polytrope above 1.64p„,i. 



energy density (^) 



pressure (^^) 



(4) 



153 
160 
180 
200 
220 
240 
260 
280 
300 
320 



2.38 
2.76 
3.94 
5.17 
6.43 
7.86 
11.27 
18.97 
30.80 
48.46 



exceeds the speed of light near the stellar center. Hence 
model stars with the APR equation of state used up to 
twice normal nuclear density and a polytrope at high den- 
sity are ruled out at the 90% confidence level. High den- 
sity polytropes that are at least moderately stiff (F > 3) 
are consistent with the 99% confidence level bound. The 
model using the APR equation of state at all densities 
is also consistent with the 99% confidence level bound, 
M/R < 0.183, but is ruled out at 90% confidence. 

Model stars satisfying the 90% confidence level com- 
pactness bound can result if the transition energy den- 
sity is chosen to be 1.64p„„ or less. For p > 1.64p„„ the 
speed of sound at the stellar center exceeds the speed of 
light for all model stars having an adiabatic index value 
large enough to yield a compactness below the 90% con- 
fidence level bound, and therefore are ruled out. For 
p < 1.64p„„ the 90% confidence level bound can be met 



without the sound speed exceeding the speed of light. 
Figure 2 shows neutron star mass vs. radius curves for 
the case where p = 1.64p„„. The radius is near 12.7 km 
for a 1.4 Mq star for the model families having F near 7, 
and the sound speed is subluminal throughout the star 
provided F < 7.023. For F = 7.023 the sound speed 
just reaches the speed of light at the center of the star. 
The central energy density and central pressure of the 1.4 
Mq, 12.7 km star are 317 MeV/fm^ and 45.1 MeV/fm^, 
respectively, when p = 1.64p„„ and F = 7.023. 

Table | lists particular values of the pressure as a func- 
tion of the energy density for the p = 1.64p„„, F = 7.023 
equation of state. Model stars that have M — 1.4M0, 
R — 12.7 km, and subluminal sound speeds are also ob- 
tained for p < l.GApnn- For p = 1.60/9„„, F = 6.25 in 
order to make the 1.4 Mq star have a radius of 12.7 
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FIG. 3: Comparison of four equations of state, three that 
use the APR equation of state at low density and a high 
density polytrope (labeled with the value of the transition 
energy density), and the fourth that uses the APR equation 
of state at all densities (labeled APR). The adiabatic index of 
each polytrope is chosen to yield a 1.4 Mq, 12.7 km neutron 
star: P = 7.023 for p = 1.64p„„, P = 6.25 for p = 1.60p„„, 
and P = 5.55 for p = 1.55/9„„ . The equations of state utilizing 
a high density polytrope differ little from each other but are 
significantly stiffer than the APR equation of state at high 
density. 



km. The 1.4 Mq, 12.7 km star has T = 5.55 when 
p = 1.55/9„„. Figure |^ compares the equations of state 
for four cases: p = 1.64p„„ with F = 7.023, p — 1.60p„„ 
with F = 6.25, p = 1.55p„„ with F = 5.55, and the APR 
equation of state used at all densities. The three equa- 
tions of state that use a high density polytrope differ very 
little from each other. Hence the equation of state of Ta- 
ble I can be regarded as an estimate of the minimally 
stiff equation of state of neutron star matter satisfying 
the 4U 1636-53 90% confidence level neutron star com- 
pactness bound that also maintains a subluminal sound 
speed throughout the stellar interior. A stiffer equation 
of state is implied if the 4U 1636-53 neutron star is ac- 
tually more massive than 1.4 Mq. Figure ^ also shows 
the equations of state utilizing a high density polytrope 
are significantly stiffer than the APR equation of state 
at high density. 



observations of the low mass x-ray binary 4U 1636-53. 
The validity of this equation of state rests on three main 
assumptions: the actual compactness of the 4U 1636-53 
neutron star does not exceed the 90% confidence level 
bound, the mass of the 4U 1636-53 neutron star is near 
1.4 Mq, and the APR equation of state is valid up to 
near 1.64 times normal nuclear energy density. The APR 
equation of state used at all densities is sufficient to de- 
scribe the structure of the 4U 1636-53 neutron star if 
its compactness is actually near the 99% confidence level 
bound, M/R = 0.183, rather than near the 90% con- 
fidence level bound value of M/R = 0.163. The APR 
equation of state used at all densities is also adequate to 
meet the 90% confidence level bound if the 4U 1636-53 
neutron star is 1.27Mq or less. (The APR equation of 
state is seen from Fig. ^ to yield a radius of 11.5 km for 
a mass of 1.27 Mq, making the compactness of this star 
M/R = 0.163.) If the mass of the 4U 1636-53 neutron 
star is larger than 1.4M0, the equation of state of neu- 
tron star matter must be even stiffer than the Table || 
equation of state. The APR equation of state is based 
on a large body of experimental evidence at low density, 
so it is reasonable to trust its validity up to p ~ 1.6p„„. 
Figure | shows the equation of state derived here (Table 
^ is relatively insensitive to the exact transition energy 
density value as long as p is near 1.6p„„. 

If the APR equation of state is used up to 1.55p„„ < 
p ^ 1.64p„„, a stiff polytrope, having an adiabatic in- 
dex of 5.5 F ^ 7, is needed to produce a model star 
with a compactness M/R — 0.163 without the speed of 
sound waves in the neutron star interior exceeding the 
speed of light. The adiabatic index would be small, near 
5/3 for nonrelativistic nucleons and 4/3 in the relativistic 
limit, if the pressure is due mainly to the Fermi kinetic 
energy. The adiabatic index is near 2 if the equation of 
state is dominated by static two-body interactions, near 
3 if static three-body interactions are most important, 
and in excess of 3 if repulsive momentum-dependent in- 
teractions are dominant | [lO| . The large adiabatic index 
of the equation of state at high density implied by the 
4U 1636-53 neutron star 90% confidence level compact- 
ness bound suggests the APR equation of state under- 
estimates at high density the strength of the repulsive 
momentum-dependent interactions. 



III. DISCUSSION 
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